Synthetic Topological Qubits in Conventional Bilayer Quantum Hall Systems 
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The idea of topological quantum computation is to build powerful and robust quantum comput- 
ers with certain macroscopic quantum states of matter called topologically ordered states. These 
systems have degenerate ground states that can be used as robust "topological qubits" to store 
and process quantum information. However, a topological qubit has not been realized since the 
proposed systems either require sophisticated topologically ordered states that are not available yet, 
or require complicated geometries that are too difficult to realize. In this paper, we propose a new 
experimental setup which can realize topological qubits in a simple bilayer fractional quantum Hall 
(FQH) system with proper electric gate configurations. Compared to previous works, our proposal 
is accessible with current experimental techniques and only involves well-established topological 
states. Our system can realize a large class of topological qubits, generalizing the Majorana zero 
modes studied in the recent literature to more computationally powerful possibilities. We propose 
three tunneling and interferometry experiments to detect the existence and non-local topological 
properties of the topological qubits. 
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One of the greatest challenges of modern physics is 
to harness the power of quantum mechanics in order to 
realize quantum computation. The fundamental obsta- 
cle is the difficulty of maintaining quantum coherence 
for long times. Among the various approaches, topologi- 
cal quantum computation (TQC)[3 is a proposal which 
relates the protection of quantum coherence with funda- 
mental properties of some macroscopic quantum states 
of matter, called topologically ordered statesjH- These 
are gapped states of matter with a certain ground state 
degeneracy determined solely by the spatial topology of 
the system. All degenerate ground states are indistin- 
guishable in all local properties, so that the quantum co- 
herence between them is robust to environmental noise, 
for a time that scales exponentially in system size. In 
two-dimensions, all low energy excitations of topological 
states are point-like quasi-particles; for a subset of them, 
called non-Abelian states, a set of quasi-particles at fixed 
positions can also carry a robust, topological degeneracy 
of states. 

The topological degeneracy in both ground states and 
non-Abelian quasi-particle states can be used as a "topo- 
logical qubit," to carry quantum information. Most pro- 
posals use the latter since it is considered easier to real- 
ize quantum operations by manipulating the quasipar- 
ticle positions. While the non-abelian states are not 
yet experimentally established, currently there are two 
types of candidate systems for non-Abelian states: The 
non-Abelian quantum Hall states[lj, the most promis- 
ing of which is proposed to be realized at filling v = 
5/2 ii, and topological superconductors with Majo- 
rana zero modes jsMTil. The topological superconductor 
proposals have the advantage of possibly larger energy 
gaps and well-controlled location of non-Abelian defects. 
However, the non-abelian defects that can be realized 
in these proposals are restricted to Ising anyons, which 
have limited computational power and the required 
proximity between superconductor and semiconductor or 



topological insulator adds complications to probing the 
non-Abelian anyons. 

In this paper, we propose a new experimental platform 
that realizes a wide class of non-abelian defects 
vastly generalizing Majorana zero modes, using conven- 
tional Abelian bilayer quantum Hall states. Our pro- 
posal is ins pire d b y recent progress in vertical field effect 
transistors, [20[ l21| and recent theoretical developments 
regarding extrinsic defects in FQH states. Il3 As is 
illustrated in Fig. [1] our setup only requires the simplest 
bilayer FQH states that have been observed in the lab- 
oratory with certain configurations of top and 
bottom gates. The key idea is to use the gate config- 
uration to induce inter-layer tunneling and thus build 
"staircases" that coherently connect the two layers. The 
end points of such "staircase" lines are non-Abelian de- 
fects, which we refer to as genons, although the bilayer 
FQH state itself is Abelian. The non-Abelian statistics of 
genons are determined by the parent state of the bilayer, 
which include the Ising anyons realized in topological su- 
perconductors as a special case. With a larger gap than 
known non-Abelian FQH states and well-controlled po- 
sition of non-Abelian defects, but without requiring su- 
perconductivity, our proposal combines the advantage of 
both types of proposed systems in the literature. 

In what follows, we explain the proposed experimental 
setup and discuss two novel experiments that can be used 
to detect the topological qubit. 
Experimental Setup 

Our experimental setup (Fig. [Ij consists of a double- 
layer quantum well system, with line junctions in each 
layer that are offset in the lateral direction. One way 
to create the line junctions is with top and bottom 
gates; the gates are used to deplete the layer closest 
to them, and to increase the interlayer tunneling over a 
small region in their vicinity. Another way to create the 
line junctions, which has been employed successfully in 
single-layer systems, [27| is through a physical obstruction 
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FIG. 1. (a) Cross-section of proposed device. Top and bot- 
tom gates partially deplete the electron fluids, and the fringe 
fields cause the wave function at the edges to extend further 
vertically, enhancing the tunneling. Filled circles and crosses 
indicate the gapless chiral edge states, moving out of or into 
the page, (b) 3D view. Interlayer backscattering between 
edge states smoothly connects the layers, allowing for coher- 
ent quasiparticle propagation between layers. The red and 
blue balls connected by a dashed line illustrates an inter-layer 
path of the quasiparticle. 



within each layer, such as an undoped A^Gai-ajAs/AsAs 
barrier. When the layers are well-separated, each layer 
can individually form its own single-layer FQH state, 
such as a filling fraction = 1/3 Laughlin FQH state. 
With larger interlayer interactions, more exotic inter- 
layer correlated states are also possible. Here we consider 
the simplest case where each layer individually forms a 
V = \/m Laughlin FQH state. Double layer 1/3-Laughlin 
states have routinely been realized experimentally, with 
gaps on the order of several K p2| - [26| . 

At the boundary of the line junctions, there are gap- 
less chiral edge states. With the offset between top and 
bottom gates, and in a suitable parameter region that 
will be analyzed below, it is possible to induce backscat- 
tering between counter-propagating edge states from dif- 
ferent layers. The back-scattering can open a gap at the 
boundary and effectively build a "staircase" that coher- 
ently connects the two layers (Fig. [T] (b) ). Such stair- 
cases are the essential elements of our proposal, which 
enable the realization of non-planar geometries. As will 
be discussed below, the end points of a staircase line junc- 
tion are non-Abelian topological defects, a pair of which 
carries topological degeneracy and thus acts as a topo- 
logical qubit. In the following we will refer to such a 
staircase as a twisted line junction (TLJ). 

Conditions for building TLJ 

The Hamiltonian of the edge states around the line 
junctions is described by the chiral Luttinger liquid the- 



ory [2i|, H = Jdx{Ho + 5Ht), with 

Ho = rriTTVo ^ naiX'^/jUpj , (1) 

a,f},I,J 

SUt = cos(m(0i7 + ^Rj) + Oij). (2) 

/,,/ 

T-Lq and 5T-Lt are the Hamiltonian densities of the ki- 
netic/interaction terms and the backscattering terms, 
respectively. a,/3 = L,R and /, J = 1,2 denote 
the chirality and the layer index of the edge states, 
respectively. The densities tIq/ = ^dxcfiai, with 
(j)ai the chiral boson fields satisfying the commutation 
relations [(I>r/li{x), (jiR/Ljiy)] = ±i^5ijsgn{x - y). 
Tjj = iTjjle"^^''' are back-scattering amplitudes between 
counter-propagating edge states LI and RJ of the two 
layers (see Fig. [T]). The velocity of the edge modes, vq, is 
assumed for simplicity to be the same for all of the edge 
states. The diagonal entries of A, Xff are normalized to 1, 
while the off-diagonal elements of the symmetric matrix 
A parameterize the density-density interaction between 
the different edge channels. 

To realize the TLJ, the following conditions are re- 
quired: 

1) The line junction of one layer must not destroy the 
FQH state directly beneath ( or above ) it. In the setup of 
Fig. [TJ this means the gates must deplete the mobile elec- 
tron density in the layer closest to them, while keeping 
the density in the further layer such that it remains in its 
FQH plateau. Based on previous experimental studies, 
such a situation is easily possible 2^ 2^ . 

2) Ti2 ^ Til, T22^ T21, i.e. the dominate back- 
scattering term should be T12. As long as the width 
of the line junction I is significantly larger than the inter- 
layer distance d, this condition can be easily satisfied. If 
the system is cleanj^ the intra-layer tunneling has a phase 
factor Tii cx e*'^/'^ (with Ib = yJhcjeB the magnetic 
length) due to the momentum mismatch, while the verti- 
cal tunneling has no momentum mismatch. This further 
suppresses other back-scattering terms. Furthermore, 
near the gates where the density is not depleted, the wave 
functions will extend further near the other layer, be- 
cause the fringe fields of the gate raise the potential well 
that vertically confines the 2DEG. This exponentially in- 
creases T12, as has been observed experimentally [20i [2l| . 

3) The back- scattering T12 should generate an energy 
gap. In the absence of interactions between edge states, 
the tunnelings Tjj are irrelevant perturbations in the 
renormalization group (RG) sense. In this large 
tunneling T12 > Hvq/Ib is required to open an energy gap 
and lead to the required TLJ. However, it is possible for 
the screened Coulomb interaction between edge states to 
make T12 a relevant perturbation in the RG sense, so that 
any small T12 will open a gap and allow coherent quasi- 
particle tunneling across the junction. The Coulomb in- 

IJ ' 



teraction gives an estimate of Ajj : rmrhvoX"'^ /Ib ^ U\ 
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where t/"^ is the Coulomb energy between edges. When 
the gates are wide {I d), the only appreciable inter- 
action will be between the two edges on top of each 
other, which, ignoring screening due to the gates, is es- 
timated by rmThvoXf-2 /Ib ~ e^/erf. In this case, the 
scaling dimension of the electron back-scattering term is 



The tunneling is relevant in the 



RG sense if Aj < 2 in the clean limit, or At < 3/2 
for disorder induced backscattering (assuming Gaussian- 
correlated disorder). For realistic parameters in GaAs 
samples, vq ^ lO^m/s, di-^ ^ IOOtito, Ib ^ 20nm, 
e ~ ISeo, we see that the repulsive interaction can easily 
make Xi-^ to be order one, so that the tunneling will be 
relevant. 

4) The edges of the line junctions should be sharp 
enough in order for the edge theory description (QJi, (0) 
to be applicable. The width of the edge region w (Fig. 
[1] (a)), in which the density of electrons gradually de- 
creases to zero, must be on the order oi Ib, so that it 
can be treated as being part of the gapless edge of the 
bulk FQH state, w is set by the fringe fields of the gates, 
and is on the order of the vertical distance h between the 
2DEG and the gate, so it is preferable to minimize h. For 
this reason, in addition to GaAs double quantum wells, 
two layers of graphene, separated by a thin insulating 
barrier may be an advanta geo us system, since h can be 
made dramatically smaller '21|. 

Now that we have demonstrated the feasibility of build- 
ing the TLJ, we will discuss the realization of topological 
qubits and experimental predictions. 

Topological degeneracy: Geometrical picture 

A I' — 1/m Laughlin FQH state on a genus g sur- 
face has a ground state degeneracy [3l|. This can 
be understood on a torus to be a consequence of the 
two independent non-contractible loops, a and 6, which 
intersect only once, combined with the fractional statis- 
tics of the quasiparticles. Together, these imply the fol- 
lowing 'magnetic' algebra for quasiparticle loop opera- 
tors W{a) and W{h): W{a)W{b) = e^'''/'^W{b)W{a). 
W{a, b) measures the fractional charge inside the loop 
a, b respectively, which can be shown to commute with 
the low energy Hamiltonian in the ground state subspace. 
Therefore, the ground states form an irreducible repre- 
sentation of this algebra, which is at least m-dimensional. 
A genus g surface gives g independent copies of the mag- 
netic algebra, and therefore a topological degeneracy of 
. On a finite size system, quasiparticle-quasihole tun- 
neling around a, b mixes the different states, splitting the 
degeneracy by an energy that is at most exponentially 
small in system size. 

Now consider our experimental proposal (Fig. [1}. 
With two TLJ's (Fig. [2]), the system now has two in- 
dependent non-contractible loops a and b that intersect 
only once; the system is topologically equivalent to a 
torus with four holes cut out on the surface (Fig. [5]). 




FIG. 2. A pair of TLJs introduces two non-contractible loops, 
so the space is topologically equivalent to a torus. The edge 
states around the outer edge of the sample and around the 
gates are mapped to edge states around four holes on the 
surface of the torus. 



With n TLJ's, the system is topologically equivalent to 
a genus g = n — I surface, with 2n + 2 holes cut out on 
the surface, which will have a topological degeneracy of 
m^. The holes cut on the surface contribute chiral edge 
states with a finite size gap inversely proportional to the 
perimeter of the hole, but they do not affect the topolog- 
ical degeneracies at all. We note that the gapless edge 
states surrounding the TLJs can also be removed through 
more complicated interdigitated gate configurations (see 
Supplemental Materials). 

Since each TLJ increases the degeneracy by a factor of 
TO, each end point of the TLJ, viewed as a point defect, 
can be associated with ^/m degrees of freedom. Similar 
to the Majorana zero mode carrying V2 de grees of free- 
dom, each end point of a TLJ is a non-Abelian topologi- 
cal defect which we refer to as a genon. (T^ The resulting 
topological degeneracies can be used to realize a topo- 
logical qubit; this setup enables manipulations of such 
topological qubits by braiding the genons. The braid- 
ing process can be shown to correspond to topologically 
non-trivial surgery operations of the genus g surface. [l2| 
These mathematical operations correspond to cutting the 
space along a non-contractible loop, rotating by 27r, and 
gluing back together. 

Proposed Experiments 

1. Quantized zero-bias quasiparticle conductance 
The y/rri topological degeneracy associated with a 
genon, i.e. an end point of the TLJ, corresponds to a lo- 
calized parafermion zero mode ll|, [l^ Il4l4l7l . |32| , which 
generalizes the usual Majorana zero modes in topological 
superconductors. Zero mode operators can be defined 
at each genon, which satisfy the 'parafermion' algebra 
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< j with a proper ordering of all 
genons in the system. Physically, the operator ai is a 
quasi-particle tunneling operator between the two layers 
at the genon point, which commutes with the Hamilto- 
nian. We refer to Ref. [12] for more rigorous definition 
of the zero mode operators. 

In the case of the Majorana zero modes localized at 
the ends of a ID topological superconductor, one of its 
key signatures is a quantized zero bias conductance peak 
with conductance G = — when a normal metal lead is 



coupled to the Majorana zero mode 0, m, m. This 

is due to the perfect Andreev refiection in the pres- 



ence of the Majorana mode. In our system, the analog 
of an electron in topological superconductor is a quasi- 
particle-quasi-hole pair, i.e., a "topological exciton" with 
charge (1/m,— 1/m). Here the two numbers represent 
the charge of the exciton in the two layers. When such 
an exciton moves across the TLJ, the two layers are ex- 
changed, such that its charge becomes {—1/m, 1/m); a 
topological exciton will become its anti-particle when go- 
ing across a TLJ, just like an electron being transformed 
into a hole in a superconductor through an Andreev pro- 
cess. The topological exciton has a perfect "Andreev" re- 
flection at presence of the parafermion zero mode, which 
can be probed experimentally (Fig. [3] (a) ). By separately 
contacting the two layers, the edge states can carry dif- 
ferent currents, and the relative current Ir ~ Ii — I2 is 
the exciton current. If we switch off the TLJ (by turning 
off the gate voltages) , the two layers are decoupled and 
the voltage in each layer is proportional to its current. 
Therefore dIr/dVr = ^x- With the TLJ the two layers 
are exchanged when the edge states meet the TLJ, so that 
the two voltages are also reversed, and dIr/dVr — 
The cases with and without TLJ correspond to perfect 
Andreev and perfect normal scattering of the exciton. 

To verify that the zero bias peak conductance observed 
from the perfect "Andreev" reflection experiment is re- 
ally from a local parafermion zero mode, we can fur- 
ther consider the experimental setting in Fig. [3] (b), 
which is an analog of using a scanning tunneling mi- 
croscope (STM) to probe the Majorana zero mode. By 
contacting source and drain to the two different lay- 
ers and bringing the edge states close to the TLJ, the 
edge states carry the exciton current 1^, and play the 
role of an STM tip. It is well-known that the electron 
tunneling from an STM tip leads to a differential con- 
ductance proportional to the local electron density of 
states. By analogy, the tunneling conductance in our sys- 
tem measures the local density of states of the exciton, 
related to G(i__i)(a;, — {V^-^ _^-^{x,uj)V(^i,^i){x,—uj)). 

Vi-i{x,t) = e*'^-f'i(=^.*)-*'^i'2(^'*) annihilates an exciton in 
the upper edge of the TLJ, or equivalently, tunnels a 1/m 
quasiparticle from layer 2 to layer 1. We find (7(1,-1) 
has the following exponential behavior: G'(i,_i) (x, w) = 
f{x)e~^^^S{u;), where ^ ~ hv/g is the correlation length 
of the gapped region (with g the gap opened at the TLJ), 
and f{x) is a power-law function of x. The exponential 
localization is remarkable, since part of the edge states 
are still gapless. To understand this, consider the edge 
theory description of the TLJ in Fig. |3] (b). In Eq. ([T|), 
([2]), we used separate chiral boson fields 4'l{R)1i 4>l{r)2 
to describe the two edges from each layer. Here, for con- 
venience we describe the whole edge of the first (second) 
layer by a single chiral boson 4'i(2)(x). Therefore 0i,2(a;) 
each lives on a circle with perimeter 2L (with L the length 
of the TLJ). The edge of each layer is a chiral Luttinger 
liquid, and it is convenient to group the two chiral Lut- 
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FIG. 3. Experimental setups to measure zero bias interlayer 
quasiparticle conductance, (a) The "Andreev reflection" pro- 
cess (see text) in which the exciton current I\ — 12 is reversed 
across TLJ. (b) Use the edge states of the sample to mea- 
sure the local quasiparticle density of states, which detect the 
parafermion zero modes at the ends of the TLJ.(c) Edge the- 
ory picture of TLJ in (b). The coordinate x increase from —L 
to 4-1/ following the direction indicated by orange, yellow, red 
and blue dots. The second layer (blue) is flipped in the right 
panel. Dotted lines between counterpropagating edge states 
indicate backscattering, generating an energy gap. 

finger liquids into a non-chiral Luttinger liquid. For this 
purpose we define the spatial coordinate x £ [—L, L] of 
the two edges with opposite chiralities (Fig. |3] (c)). x 
increases along the counter-clockwise direction for layer 
1 and clockwise direction for layer 2. With such a choice, 
(/)i and (/)2 has opposite chirality (defined with respect to 
i), and the vertical back-scattering induced at the TLJ 
occurs locally between (f>i{x) and 4'2{x). The Luttinger 
liquid Lagrangian can be written as 

L=j di[{dt^f +v{d5:ipf + ge{-i)cos{jnip)\. (3) 

with ^ = ipi -\- 4>2 aud where 9{x) is a step function. The 
tunneling induces a gap in the region x £ [— L, 0], and the 
region x € [0, L] remains gapless. In this representation, 
the m topological degenerate ground states are correctly 
given by the m minima of the cos{mip) term, with ip 
pinned to the values 27rn/m, n = 0, 1, m — 1. Impor- 
tantly, in this coordinate the exciton annihilation opera- 
tor becomes nonlocal: Vi,_i(x) = e*'^i(^)~"^2(~^'. Due to 
the non-local commutation relation of chiral boson fields, 
applying this operator to the ground state causes a do- 
main wall in if in the region ~L < x < 0. Away from the 
defects at a; = and x — ~L ^ L, such a domain wall 
costs an energy of order the gap g. Therefore, the corre- 
lations of this operator must decay exponentially, away 
from the positions of the defects. At the positions of 
the defects, this operator keeps the system in the ground 



5 



state manifold because it does not create a domain wall 
in in the gapped region. Instead, it simply yields a uni- 
tary transformation among the ground states by globally 
flipping the value of ip in the gapped region. From this 
analysis we see that the zero bias peak in the exciton 
local density of state is indeed probing the topologically 
degenerate states. 

2. Long-time current-noise cross correlation 
Since the genons are non-abelian defects, it should be 
possible to detect them through interferometry. Com- 
pared to the interferometry measurements proposed 
for detecting non-Abelian anyons in non-abelian FQH 
states 0, an advantage of our current setup is that 

the non-abelian defects (genons) are imposed externally 
rather than being dynamical quasiparticle excitations, so 
that their position and number can be easily controlled. 

Our setup (Fig. |4]) consists of two TLJs and four quan- 
tum point contacts (QPCs), created by either a top gate 
or a bottom gate. One pair of QPCs (ri^2 in Fig. al- 
lows for weak quasiparticle tunneling across the normal 
region of the fluid, between the edges of the top layer. 
The other pair (ri^2) allows for weak quasiparticle tun- 
neling across the TLJs, from one edge of the top layer 
to the other edge of the bottom layer. It should be no- 
ticed that there are only two nontrivial interference loops 
formed by the four QPCs, the first one (L in Fig. |4|) in the 
top layer including Fi and and the second one (L') 
going through both layers including Fi and F2. Other 
interference paths such as the one enclosing Fi and Fi 
does not form a closed loop since the quasiparticle can- 
not return to the same layer. The interference in these 
two loops L, L' measure the electric charge in the loops, 
which measure the topological qubit state of the genons 
1,2 and 2,3 respectively. As a key consequence of the 
non-Abelian nature of genons, the topological qubit de- 
grees of freedom are non-local, and the topological charge 
carried by genons 1, 2 does not commute with that of 2, 3. 
Therefore non-local correlation between interference pat- 
terns in these two loops are induced due to the uncer- 
tainty principle. Alternatively, the non-commutativity 
of the two quasiparticle paths L and L' can also be un- 
derstood from the geometrical picture shown in Fig. [21 
with the measured topological charge corresponding to 
W{a), W{b). 

To be more precise, the interference pattern can be 
observed by applying a voltage difference V between the 
left and right-moving edges in both layers, and mea- 
suring the currents /i, I2 (Fig. |3]). We find that 
the interlayer current noise cross-correlation 5*12 (i) = 
i({/i(t),/2(0)}) - (/i(t))(/2(0)) has a dramatic feature: 



5'i2(|i| > 1/T) ^ constant ^ 0. 



(4) 



Such a long time noise correlation, at finite temperatures, 
is a striking direct consequence of the non-Abelian nature 
of the genons, which has never been proposed or observed 
before in an Abelian state. Intuitively, this is understood 



FIG. 4. QPC interferometer setup for detecting the non- 
Abelian nature of genons. ri,2 are the two QPCs in the 
first layer, and ri,2 are the two inter-layer QPCs across the 
TLJs. L and L' indicated by thick dashed lines are the two 
interference loops. 



as follows. The current Ii (t) depends on the quasipar- 
ticle interference between the two QPCs Fi.2- When a 
quasiparticle leaves the top layer through the QPC Fi, its 
fractional statistics cause a phase shift in the interference 
of loop L and therefore changes the value of Ii{t). Such 
a phase shift is permanent until the next quasi-particle 
tunneling occurs at Fi. The converse is also true, that the 
quasiparticle tunneling at F2 causes a permanent phase 
shift in l2{t)- Therefore, the fiuctuation of the currents 
Ii{t) and 12(0) have a long-time correlation for both t > 
and t < 0. It should be noticed that only having three 
QPCs Fi,ri,F2 [? ] can already induce a long-time 
correlation of 5i2, but only for t > 0. 

We present in the Supplemental Materials an explicit 
calculation of 6*12 based on the edge theory; the result, for 
1^1 ^ 1/T and to lowest order in the tunneling strength 
at the QPCs, takes the form: 



Suit) ^T^-^[\Tlr,t2\0{~t)s\l^ + \rlr,r2\e{t)s\t 



T,T2r,T2\{s[feit) + s\ye{~t))]. 



(4)/ 



(5) 



The dimensionless functions S'12 , whose explicit form 
is given in the Supplemental Materials, depend on the 
quasiparticle statistics (they are non-zero only when m =/= 
1) and decay with the dimensionless ratios HeV/mkBT, 
and ksTx/Hv, which are optimally of order one, where x 
is the distance scale between the QPCs. 

On much longer time scales, thermally excited mobile 
quasiparticles can decohere the state of the qubit, which 
will wash out the noise signature considered here. How- 
ever this time scale is expected to be exponentially large 
in 1/T, (ie t oc e^/-^), because the quasiparticle density 
is exponentially small in 1/T. 

Discussion 

A longer term goal is to implement unitary gates on 
the topological qubit by braiding genons. This can be 
realized with further sets of gates to tune interactions 
between various zero modes localized at the genons. [12. 
14 , 3, 13] The computational power of genons in generic 
Abelian states is stronger than that provided by Ma- 
jorana zero modes 14|. The setup proposed here can 
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be generalized to realize universal TQC in non-abelian 
states that by themselves would not be universal for TQC 
[l^ . Examples include two layers of the non-abelian 
phase of Kitaev's honeycomb model 12, 3^, or the odd- 
denominator Bonderson-Slingerland FQH states 
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FIG. 5. The defects (genons) can be viewed as domain walls between two different ways of generating an energy gap in 
the counterpropagating edge states (double arrows indicate interedge electron tunneling). The zero mode operator 021-1 
corresponds to the quasiparticle hopping process shown in the figure, projected into the ground state subspace of the edge 
theory. 



Supplemental Materials 
I. Parafermion Zero Modes 

As mentioned in the main text, an important property of the TLJs is that their ends locahze exotic forms of zero 
modes, called parafermion zero modes, which generalize the usual Majorana zero mode algebra. [l2j Here we will 
briefly recall the definition and physical meaning of these zero mode operators, reviewing results from [l2| . 

Let us imagine aligning all TLJs, and cutting the whole system along the line, yielding the counterpropagating edge 
states. Then, the normal regions of the fluid can be obtained by gluing the edges back together with the electron 
tunneling terms X]/(^eL/^e-R^ + -^-'^-) ^ cos(m(/)/), while the twisted line junctions can be described by gluing the 
edges back together with the twisted tunneling Vl/g^j^'effi + ^'ei2*efli + ^ cos(m0i), where we have defined 
(f>i = 4>Li + (l>Rij and = (j)Li + 4'R{i+i)Yo2- Here, we have included also the tunneling term vE'gj;^2^e^Ji + which 
is not present in the experimental setup of Fig. [T] 

The ends of the TLJ can therefore be considered as domain walls between these two different ways of gapping the 
counterpropagating edge states: 

,nj _ /E/COs(to07) if xeAi 
I LiCos(m(/)/) if x e Bi 

where Ai indicates the normal regions, while Bi are the regions associated with the TLJs. 



As shown in 12|, we can define quasiparticle tunneling operators (see Fig. [5]) Q!2i-i = g^^'I'^i^Ai) ^ «0i(a:B. )^ 
ct2i = e^^'^^^'^i^ e~'^^^^^^^^^+'^\ which are zero modes of the edge Hamiltonian, [ai, Hedge] — 0, and which satisfy a 
Zm 'parafermion' algebra: ataj = ajaiC^^^^"^ , and a™ = 1. 

Other ways of creating domain walls between two different wa ys of g enerating an energy gap in a single pair of 
counterpropagating chiral Luttinger liquids has been explored in [14l4l6l |. where similar parafermion zero modes are 
obtained. 



II. Interdigitated gating 

In the experimental setup described in the main text, there are gapless edge states surrounding the TLJs. These 
edge states do not affect the topological degeneracies, the exponential localization of the zero modes, or the possibility 
of braiding the genons by tuning the interactions between the parafermion zero modes. Nevertheless, it might be 
desirable to remove these gapless edge states, while still preserving the topological degeneracy. This can be done 
simply by alternating the offset in space (see Fig. [6]). The tunneling term in the edge states can be modelled as 
follows: 

5T-Ltunn = a{x) cos(m(^i) + /3(x) cos{m^2)- (7) 

While a{x) and fi{x) are never both non-zero simultaneously, at long wavelengths only the average value of a and P 
is important. Therefore, both edge states can be gapped in the twisted way on average in this setup. 
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FIG. 6. Alternating the offset of the top and bottom gate will allow us to generate an energy gap for all edge states. At long 
wavelengths, both twisted tunneling terms, J2i cos(m(^/) are simultaneously present. 



III. Interlayer Current Noise Cross-Correlation 



A. Summary and Results 



Here, will present the derivation of the current noise cross-correlation for the interferometer summarized in the 
main text. We consider a tunneling Hamiltonian with four QPCs (4ll. |4^: 

2 

i=l 

H,{t)=i,d+{x.,)+i*d-{xi), (8) 

where 

v,=T,e-'^'\ i,=f,e-'^-K (9) 

Ti,Vi are the tunneling matrix elements for the different QPCs, and uji — e*Vi, where Vi are the potential differences 
across the QPCs and e* — e/m is the elementary fractional charge of the quasiparticles. We note that in the setup of 
Fig. m to lowest order in the tunneling matrix elements, the voltage across each QPC is Vi = V . For what follows, it 
will be useful to define: 

a;i2 = a;i-X2, Ji2 = ii-i2- (10) 

The operator is the quasiparticle tunneling operator across QPCs enclosing the normal regions of the fluid: 

O±(x,) = ^±e±^(*--*i«)(x0, (11) 

where Hi are Klein factors that ensure that the quasiparticle tunneling operators at different points in space commute 
with each other, as required. The operators are quasiparticle tunneling operators across the twisted line junctions 
(TLJs): 

a±(i,0 = ^;fe±^('^--*^«)(iO, (12) 

where nf are also Klein factors. The commutation relations of and can be determined from the boson 
commutation relations: 

TT 

[4'Li{x),<j)Li{x')] = i — sgn{x - a;') 
m 

TT 

[0i?j(a;),0j?j(a;')] ^ -i — sgn{x - x'), (13) 
m 

from which we determine that [k^, = 0, [ki, k|] = 0, and 

K^R^ = K^K^e^'"^^ , (14) 

where 9^^ = rsnsgn{xi — Xj)/m. 
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In the interferometer discussed in the main text, which probes the non-abehan character of the genons, we require 
xi < ii < a;2 < £2- However it will also be helpful to compare this with the case cci < a;2 < ii < i2, which does not 
probe the non-commutativity of the two non-commuting, non-contractible quasiparticle paths. 

Now, we let /bi and 7^2 be the current that is backscattered through the normal regions and the TLJs, respectively: 

h.j{t)^^^^-l[QRJ,Htunit)], (15) 

where Qbi = J dxdx4>Ri is the total charge on the right edge in the /th layer. Thus: 

2 .2 .2 

Ibi{t) = - V[r,e-*^*^^*O+(x0-i/.c.] = -y\n,0+{x,)-v*0-{x,)\ ^ -Y^Y] €,v1:0'^{x{) 
m ^ — ' m ^ — ' m ^ — ' ^ — 

i—\ i— 1 i—\ ei—± 

2 .2 .2 

h2{t) = A V[f,e— *^=*d+(x,)-i^-c.] ^-y\i,d+{x,)-v*d-{x,)\ = - V V e,v^d'^{xi) (16) 

TO ^-^ TO TO 

i=l i=l i=l ei=± 

In the setup of Fig. 01 we find that: 

1 e2 

h = -^V + Ibi, 
TO n 

-^V + h2+Ig2, (17) 

TO, h 

where Ii and I2 are the currents flowing through the device in the first and second layer (see Fig. |4]). In order to 
understand 7^2, observe that current can be backscattered through the TLJs in two different ways. The first way 
is through quasiparticle tunneling across the TLJ. The second way is if the outer edge that follows the constriction 
of the QPC overlaps the gapless edge states that flow around the TLJ. In this case, there can be interlayer electron 
tunneling between the top and bottom outer edges of the sample, through the edge states surrounding the TLJ. This 
latter contribution is what we call Ig2. Therefore: 

= i(/bl(i)(/fc2(0) + 1,2(0))) - (/l(t))((/2(0) +7,2(0))) (18) 

We will be interested in the limit \t\ ^ l/T, in which case it will become clear that 7,2, which involves only electron 
tunneling, will not have long-time cross-correlations between the quasiparticle tunneling current Ibi- 
Therefore, we will compute the following interlayer noise cross-correlation: 

Si2it) = i({7w(i),7fc2(0)}) - (7bi(0) (7,2(0)). (19) 

Note that the expectation values are given by 

(Ait)) = (rcA(Oe-'/o^-"(-)'^-)o, (20) 

where the (■ ■ • )o indicates an expectation value in the unperturbed theory. C indicates the Keldysh contour, and Tq 
is time-ordering in the Keldsyh contour. In our calculation, we will drop the subscript and assume all expectation 
values are in the unperturbed theory. 

Let us first present the result of the calculation. We find that to lowest oder in the QPC tunneling matrix elements, 

Si2{\t\ » i/T) ^ (^^)'^;" \ \r,\^\rir2\9{-t)sg\cj2,ioi,xi2,^i2) + \ri\^\rir2\0{t)s['^{oji,u;2,xi2,vi2)+ 

TO-^ 

|rir2f if 21(5^2^ (tlJi, ^2, X12, Xi2,Vl2,Vl2)0{t) + 5^2^(0:2, Wi, ii2, Xi2, ipi2, (pi2)0{-t))\. (21) 



where we have defined the dimensionless quantities 



h LOl _ _ h UJ2 

IT Kb 1 ttkb -I 

- _ T^ks _ 7rfci3 

Xi2 = -^Txi2, Xi2 = -^—Txi2 (22) 
nv nv 
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and tpi2 and (pi2 are the phase differences between the tunneling matrix elements of the QPCs: rir2 = |rir2|e''^i^, 
and f = |f if 216*"^^^. Note that 5^2^ depend on only dimensionless quantities. 

S'^2'(w2,wi,a;,(^i2) =4c,.((:j2,2a;)[(cos(27r/m) - 1) cos(^i2/_(wi, 0) - sin (^12 sin(27r/m)/+((Iii, 0)], 

{C)i,u)2,x, y, <^i2, <y5i2) =8cr ((Iji , 2a;) cos if 12 [2{cos{TT/m) - 1) cos(<^i2)cr((:j2, 2j/) 

- sin(7r/TO) sin(¥3i2 )/+ {ui2 ,y)], (23) 

The functions f±{uj,x) are defined as 

\sinh |i + sinh |t — i| / 

while c,-(a, h) is defined as 

c.(a,.) . -2sin(.M) Im ..^^ ^ ^'"} (25) 

As a check, observe that when m — \ and to — > oo, the two non-contractible paths for quasiparticle propagation 
commute, and therefore this noise cross-correlation should vanish at long times, which it indeed does. 

In Fig. [71 we plot the functions f±[a,h) and —Cr{a,h)/ sin(7r/m). In FigEl we plot S'Jj (a, 6, a, 6, 0, 0), for i — 1,2; 
that is, for the case where the voltage differences between all QPCs is the same, Xi2 = Xi2, and for zero phase differences 
(pi2 and (pi2. We see that when a and b are of order one, reach their maximum value, and decay quickly as a and 
b increase. From this, we see that the noise cross-correlation will require the dimensionless parameters w^, Xi2-, and 
X12 to all be of order one. If we assume T « 50 — IOOto^C, and v w 10^to/s, we see that this requires X12 and X12 be 
on the order of several fj,m, which is reasonable for realistic device dimensions. 



B. Klein Factors 



Here, we will determine the Keldysh-contour ordered correlation functions of the Klein factors, which are required 
for the calculation. To do this, first we observe that the two point function for Klein factors that commute with each 
other should be one. Next, we observe that the correlation functions of several Klein factors should factor into two- 
point functions. In order to obtain the two-point function of non-commuting Klein factors, it is helpful to represent 
the Klein factors in terms of bosonic operators 43| : 



(26) 



where 9i, 
fact that 



6i are bosonic operators and for what follows we will leave the normal ordering 



: implicit. Using the 



:: : e^^'^l for cases where [A, B] is a c-number, we find that the commutation relation of 




FIG. 7. Density plots of (a) -j^^fe^, (b) /-(a, 6), and (c) f+{a,b), for m = 3. 
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(a) 




43.9 



(b) 



FIG. 8. Density plots of (a) s'^^a, b, a, b, 0,0) and (b) -S[2' {a,b,a,b,0,0), for m = 3. Both are zero when a = b = 0, reach a 
maximum for a and b of order one, and decay to zero with increasing a and b. 



the Klein factors is reproduced when 

[9i,0j] = -iiTsgn{xi - xj), [Oi,d]] = = 

Since 9i and 9j are conjugates, we can form a harmonic oscillator algebra: 



1 



+ it 



1 



2n 



such that 



[at] , a-j \^ ~sgn{x^- Xj), (a^ a ^ ) = . 



(27) 



(28) 



(29) 



We suppose that the Hamiltonian for these bosonic fields is zero, as they have no dynamics. From this we conclude: 

~ ~ - - 2,71 

{0^{t)e,{o)) = {e^{t)9,{o)) = tt/2, {e^{t)0j{o)) = -{ej{t)e,{o)} = -—sgn{x, ~ xj). (so) 



From this, we conclude that the Keldysh contour ordered correlation function is: 



where 



Xmv2W = — — S9n[t) — . 



Therefore: 

(rcK-(t''^)K^(0''^)) = e^<^^^Xviv2W''9n{x,-Xj) 

For later reference, we will need the following correlation function: 



(31) 



(32) 



(33) 



(34) 

To calculate this, we use the fact that (Tce"^^ • • • e"^") — e^'<j^"^'^"^'"*^\ where in the Klein factor case, we only include 
correlation functions between non-commuting variables. Therefore: 



(TcKq {t^)K'^'^^ {t'^^)k'^^ (0^'')k^'^^(0''^)) = e^^'^^'^2{Xv-vW0o,ia,2-X,l,12(*)'^''lfi2-Xvi-vW'^ftlc.2+Xvi12W^f3lft2) 



e~''^isr'^[''(''°l=2-'*/9ic2) + '?2(Sai/32-S3i/32)l t>0 
e^T^'^[''(*°l'32^«<»l°2) + '?l('*/3lC2-«/3l^2)l ^ < 



(35) 
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where Sap = sgnixa — xp). When xi < 2:2 < aii < X2, Sap = —1 for all a, /?, and therefore 

(rc<(t'')^^;Hi''^)^a^.(o-'')^^;^(o''^)) = 1. (36) 

When xi < ii < a;2 < 2:2, which is the case of interest for our proposed interferometer, we instead have 

C. Boson correlation functions 

We will need the following correlation functions: 

G,,,,-(x-x^^-^0^ (rce-'^(-'*'')e-"'^(^''*'''')) = e<^^^(^'*')^("''*''^ (38) 
At finite temperatures, this is (setting the velocities of the left- and right- movers to unity): 

G ut)=( (39) 

where 

Xmn2W = — 2 — ^^^^ > 2 — ^ ^ 

Taking the limit (5 — > 0+, we can replace the above by 

U^,^,[x,t)-e {smh[7rT\t + x\]smh[7rT\t-x\]J ^ ' 

The chiral correlator, involving only the left or right-moving modes (and re-inserting the velocities Wi/i? for the 
left/right movers), is: 

i^,,,A^,t) - [Ice )-e [smh[nT\\vL/n\t±x\]) ^^'> 

D. Details of calculation 

Here we describe in detail the calculation of Si2{t), to lowest order in the tunneling matrix elements. We have 
(recah eq. ^ ): 



\{{hl{t),Ib2{0)}) = ^ ^(-*)' / dTidr2(Tc/bl(f' )^rt.n(n)/b2(0-")i?t„„(T2)) 

1 f°° 

= -0 / Vimdrrdrr{Tchi{nH{Tr)h2{0-'')H{rr)), (43) 



'72 =± 

Consider the correlation function: 



ai, 02,01,02 «ai«^iea2<^|32 

(TcO^- (x„, , fr)0'^^ (a;^^ , )o^^2 , (j^^ , )) (44) 

The above correlation function vanishes unless = ^e^Ji and — Therefore: 

(Tc/w(i'')i?(r?^)/b2(0-'')7?(r2^^)) = - -1 ^ E<V.''<2«r^i^2(Tc<(t'')^^;H^r)^a^2(0"'')^r(^2^)) 

0-1,012, Pi, 02 ^1^2 

^J|^g-iei</'lH(^c,i,t'Og»ei0iJ?(a:/Ji.'ri''^)^ ^J|^g-*«2</'2H(5Q2'0~'')e"2</'2H(i32''^™ 
(Tce"^*^^-^*-^"! ,t'')g-«l0ii,(x^j ,r;"-)gie2</>ll,(a:c,2,O'')g-ie20lL(/92,r2'2)^ j-^j^-j 
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Recall that in our notation, the superscript implies complex conjugation, not inverse. 



ei,e2— ± ai,(3i,a2,P2 

^J'pgiei0i-i(a:ci,t'')g-i«i0i-i(a;3i,Ti''i)gj£20ii(xa2'O~'')e~*^20i,i(5/32''^2^)^ (46) 

This gives: 

{Tchi{f')Hir'nh2{0-')H{Tr))=-^ E E £ie2r^\r-T^^^r-^e-^-^(*--^)+-=-=(--=) 

(rc<\ (rr)«^^2 (o-'')«^;^ (rr)) 

^-'vm (^"2 ~ ^/92 J ~'''2)G_'^^^ {Xa2 ~ ^02 ) "'''2) 

/ G^^r"" (^"1 - '^/32 : ^ - T-2)G^;-^ (a;/3i - , ri) ' 



(47) 



The Green's functions G^^, (x, t) decay exponentially in t; therefore, the integral will only be dominated by regions 
where ti « t and T2 « 0. Furthermore, we are interested in the long time behavior, where t is large. In this limit, the 
ratio in parantheses above will go to one: 

/ G^pn'^ (a^ai -X^^t- T2)G^;C^^ {xp^ - , n) \ _^ ^ ^^^^ 
\ G^l^^"^ (Xaj — 1 i)Gniri2^ {xp^ — X , Ti — T2) J 

when t « n — >■ 00 and \t\ ^ \t2\, t ^ \xa — Xjs]. Therefore: 



ei,e2=± ai,/3i,a2,P2 



G%«^^^- - i^, , -r2)G;;,f (a:„, - x^, , t - n) 



(49) 

Therefore, so far we have, for large t, 

';i';i,'72ei,e2=± ai,/3i,a2,/32 

[(Tc< {t^)K-pl' {t^^)K% (0-");^-= (0"=)) - (Tc< (t")/.^;^ (t"^ )) (TcK- (0-"))^-= (0"=))] 

/OO /'OO 
gie2a,2(-r2)G--;ji. (j^^ _ i^^, / drie-i-i(--^)G;;;,f''- (x«, - x;3, , i - n) (50) 

-00 J —00 

In the case where xi < X2 < xi < X2, using the formula for the Klein factor expression found earlier, we find that 
the difference of Klein factor correlation functions vanishes, so that 6*12 (t) — for \t\ ^ 1/T. 
Now let us focus on the case xi < xi < 0^2 < 52- In this case, for t > 0, 

> 0) E E eie2r^^^r-^f^=^f-ne-^'^("^-''^^^[''^~-^+''=^^>=^^J - 1] 

T7,'?i,T?2ei,€2=± ai,/3i,a2,/32 



/•OO 

mm 



/OO /'OO 
ei.20.2(-r2)c,.2H;^.iz,(j^^ - -r2)(ir2 / dne"--^^--^)^;;^,^ (x„, - - n) (51) 
-00 J —00 



In what follows, let us assume that the edge velocities are the same on the left and right edges, and in the two 
layers, which we set to unity: vil = vm = V2r = 1- 
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Now let us define: 



\smh[TTT\t + x\]smh[TrT\t-x\]J ^ ' 



Note that this satisfies 

A-^{x)=Az-4x), A^^^ix) = A-^^^{-x), A-^{x)+A-_{x)=A-_{x)+A%{x). (53) 

Then we can write the noise as a simple expression: 
1 1 
2 



J),?)l,I)2Cl,e2 = ± Ql,^l,Q2,/32 



and 

^i2(i<o) = i-^ E E eie2^i^2|r„,r^,f„,f^je"^^-''i+"^'^-^^ 

'),')l.»;2ei,e2 = ± ai,^i,a2,)92 

[e-.^(o2-/52).[-.^„,2+m^.,.] _ l]4^,;-i(x„,0jA!_^-^^(i„,^,) (55) 

Now let us try to make more explicit the above expressions. We start with the t > case. Summing over a\ and 
/3i gives: 

5i2(t > 0) =i^|rir2| ^ ^ eie2|fa2f;3,|e*^i'^"+^^^'^°^^Ke"^''™^''''°^''*"''^"^'' -1) 

mmA'^:;^^ (ii2)(^^i,? (2:12) - A-;^"^ (0:12)) (56) 
For 5*12 (t < 0), we perform the sum over a2, 1^2 to get: 

^i2(i<0) = i^ ^ ^ eie2mr?2|r„,r^,fif2|e"^^-''i+"^^- 

'),'Jl,»72ei,e2=± ai.^i 

[e.^.[-.^.,2+..^.,2i _ i]A--(^o,ft)(Ai^,%^,(5;i2) - ^:',C(ii2)) 

(57) 

This can be written as the sum of two terms: 

^12 (i > 0) = 51^2 {t>0) + Sf^ it>0), (58) 

where 5*^2 (^ > 0) contains the terms where a2 ^ /32 in the above sum. Sf2{t > 0) contains terms in the sum ([55)1 
where a2 = 132- Similarly, we write 

^12 (i < 0) = 51^2 (i < 0) + 51^2 (i < 0) , (59) 

where S\2{t < 0) contains terms in the sum ([57]) where ai ^ while S'i2(t < 0) contains terms in the sum where 
ai = After simplifying, we find 

/j^rji\A/m-2 _ _ 

S'i2(* > 0) = 5 |rir2fif2|S'i2 ('^1,^2, 7rTxi2,7rrii2,<^i2, 9312), 

^?2(* > 0) = \V^\^\T^T2\SY^{UJI,UJ2.^TX^2.V12), 



m — 2 



5'i2(^ < 0) = „ |rir2fif2|S'^2\w2,wi,7rTii2,7rra;i2,'^i2,<^i2), 

^?2(i < 0) - ^—^2 — |rip|rir2|^l2H<^2,<:si,7rTii2,(^i2), (60) 



where are presented in (f^. 



